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1 Introduction 

Exiting 50 years of foundation of Precision EW Physics were perfectly re- 
viewed in several brilliant talks at this Symposium £Q. I hope, our commu- 
nity would share my opinion that after 12 years of excellent work of LEP we 
have full rights to say that a new scientific discipline has been born, Precision 
High- Energy Physics, PHEP, consisting of: 

— experimental measurements themselves at per mill precision level; 

— supporting theoretical calculations with even better precision. 

Me and my colleagues were involved into this glorious period with ZFITTER 
project. Actually, ZFITTER project has been started 25 years ago, in 1976, 
when we wrote our papers on on-mass-shell renormalization scheme. For- 
mally, first ZFITTER was born within CERN workshop "Z Physics at LEP 1" 
in 1989. As an important intermediate step I would like to mention CERN 
workshop 1994-1995, "Precision calculations for Z resonance", which I had 
pleasure to run and to have Alberto among our contributors. In 1999 we 
have published a long write-up of ZFITTER v6.21, which recently appeared 
in CPC [2 . When I say "12 Years of Precision Calculations for LEP", I 
mean that during all lifetime of LEP, 1989—2000, ZFITTER was upgraded 
and supported in ADLO, SLC and LEPEWWG. We would like to see the 
end of ZFITTER project in the Y2K and, therefore, a very natural question 
arises: What's next? 

Meantime, me and Giampiero Passarino in 1997-1998 wrote a monstrous 
book: "The Standard Model in the Making", which appeared in 1999 in 
Oxford University Press It has some value for this talk. 

ZFITTER uses an one-loop core and general environment based on our own 
formulae, and incorporates all the world results for higher order QED, QCD 
and EW Radiative Corrections (EWRC). As an input it uses the so-called 
LEP1, IPS - Input Parameter Set, i.e. 5 parameters: 

Aa£\M 2 z ), a s {M 2 z ), m t , M z , M H . (1) 

M z is measured very precise at LEP1 and for the first three parameters a 
rich information is available from the other measurements. Therefore, we are 
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approaching one-parameter fit with M H being the only parameter. Results 
of such a fit are being presented in famous Blue Band figure, derived with 
the aid of TOPAZO & ZFITTER codes. 

ZFITTER incorporates all known in the literature QED RC up to O ((aL) 3 ) 
(where L = Ins jm\ — 1 = 23 at s = M 2 and, therefore, effective QED 
coupling is quite large: aL = 0.169); QCD RC up to O (a 3 ); and EWRC up 
to O (a 2 ), however, for the latter only leading and subleading RCs are known 
(Degrassi, Gambino, Sirlin et ai, 1996-1997). 

In several cases perturbative calculations are saturated like it takes place 
for leptonic contribution to the running QED coupling, Aa;(s), which is 
presently known up to three loops (2-loops Kallen, 1955; 3-loops Steinhauser, 
1998): 

A ai = 314.97637-10- 4 = [314.18942 1 _ loop +0.77616 2 _i oop +0.01079 3 -i O op] -KT 4 . 

(2) 

The last known term is small enough to be used as an estimator of theoretical 
uncertainty. Typical scales of the problem: LEP1,2 energies \/s = M z — 
200 GeV , M w , M z , mt , M H , all are of the same order and calculations must 
be complete and one-loop calculations were complete! For two-loops, the 
notion of mf enhanced terms was introduced. Only terms O (G F mf) and 
O (G F m 2 t M 2 z ) are known. Since likely, 100 GeV < M H < 250 GeV, popular 
expansions in M^/m 2 and rn^/M 2 have bad convergence. One may say 
that complete two-loop EWRC were welcome for LEP1 but our 
community did not deliver them! 

In this talk I will tell about three directions of possible future work in 
which we were digging in the Y2K. The first one we call: 

2 Book "heritage" 

While working on the book, me and Giampiero Passarino wrote dozens of 
book supporting form codes. Later on an idea was erupted to collect, order, 
unify and upgrade these codes up to the level of a "computer system": to 
which we gave the name CalcPHEP - - 'Calc'ulus of Precision High Energy 
Physics. 

This system is being realized at a web site: brg.jinr.ru. Presently 
we do not maintain any author list, we say that we have people who con- 
tributed to it: G. Passarino, DB, L.Kalinovskaya, P.Christova, G.Nanava 
and A.Andonov [3]. 
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After completion of R&D phase of the project, presently the following 
options are available at the site: 

1. generation and reduction to the scalar PV functions of one-loop Feyn- 
man diagrams (FD) for all SM 1 —>■ 2 decays and 2f — > 2f processes 
(in i?£ gauge, QCD included); 

2. computation of one-loop scalar form factors for decays Z(H) — > // and 
W -> /;/,: 

The system was already used for calculation of EWRC to e + e~ — > tt (as 
described in Section 3). 

In case if the project will be approved by scientific community (I mean, 
it will be decided that it is worth doing, since there is certain competition 
with FeynArts, see, for instance j^j), during next two years (2001-2002) it 
would be feasible to realize next steps: 

1. extension of computation of one- loop scalar form factors for all SM 

1 — ► 2 decays and 2f — ► 2f processes of experimental interest (REI 
criterion - Reactions of Experimental Interest. By this we mean that 
in the initial state could be only e^,^, muons, neutrinos and partons). 

2. extension of availability of generation and reduction of FD for all 2 — > 2 

processes; 

3. realization of the step: from form factors to helicity amplitudes; 

4. realization of the step: from amplitudes to realistic observables; 

5. solution of the problem of automatic generation of codes for numerical 
calculation of form factors, amplitudes and observables; 

6. creation of codes for calculation of one-loop amplitudes for SM REI 

2 — > 3 processes; 

7. user support of CalcPHEP system; creation of user-friendly environment 
on the site. 

3 New calculation for e + e~^ tt 

The process e + e~ — > tt is studied already about ten years in connection with 
experiments at future linear colliders, see for instance recent review jS]. 

Actually, it is a six-fermion process, however the cross-section a(e + e~ — > 
tt) with tops on- mass-shell is an ingredient in various approaches like DPA [7j, 
or the so-called Modified Perturbation Theory (MPT) [8 . 
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Recently we completed a new calculation of the electroweak part of the 
amplitude of e + e~ — > ti process jU] in two gauges, and unitary ones. 

There were many studies of e + e~ — > ti process in £ = 1 gauge, see for 
example [12] and [TT] . 

The purposes of this new study are: 1) to explicitly control gauge invari- 
ance in and search for gauge invariant subsets of diagrams in fully massive 
case; 2) to compare with the result in the unitary gauge as an internal cross- 
check; 3) to propose a way of realization of the step from FD to renormalized 
amplitudes within CalcPHEP project; 4) to compare with existing in the li- 
terature results; 5) to create a FORTRAN code for the IBA (Improved Born 
Approximation) cross-section dc/dt ~ \A\ — I^-bo-rjv 

+ A WEAK \ for subse- 
quent use within MPT and within "algebraic" approach, see Section 4. 

3.1 Amplitudes in L,Q,D basis 

In presence of massive fermions it is convenient to introduce the so-called 
L,Q,D basis in which the amplitude may be parameterized with six scalar 
form factors: 

AnQ e Q f 
Ay = % — - a(s) 7m® 7m ; 

A z = ^e 2 4/( 3 )/f^{ 7 , 7+ ® 7 , 7+ F ii ( S ,t) 

-4|Q 6 |4, 7m ® 7m7+ F ql (s,t)-4\Q t \s 2 w ^ + ® lfl F LQ (s,t) 
+16|Q e Q t |s^ 7/i ® 7/1 F QQ (s,t) 

-l^l+®^rn t D^F LD (s ) t)+A\Q e \s 2 w ^^®im t D^F QD (s,t)J. (3) 

where 7+ = 1 + 75, = {pi — Pt)„ an d Xz(s) is the 7 /Z propagator ratio. 
Every form factor in R^ gauge could be represented as a sum of two terms: 

Fl,q,d (s) = F^) }tD (s) + F1 d Q D (s, £) . First term corresponds to £ = 1 gauge 
and the second contains all £ dependences and vanishes for £ = 1. We checked 
the cancellation of £'s separately for six subsets of diagrams: 1) with virtual 
7 's, i.e. QED subset; 2) — Z, 0°, the so-called Z cluster, i.e. vertices with 
wave function renormalization factors, see Fig. ^ 3) — H,(f)°, i.e. H cluster; 
4) — W, (j) , i.e. W cluster plus all self-energies and WW box; 5) four Z 7 
box diagrams; 6) two ZZ box diagrams. 
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3.2 Scalar form factors of Z cluster 

As an example we discuss scalar form factors of Z cluster, Fig. 

7/ 1, 




z + 




+ 




Figure 1: Z cluster. 

They all are gauge invariant objects and four of them are UV-finite. We show 
pole contributions (1/e) explicitly and denote finite quantities with T (s): 

FT 00 = H Z («) , (s) = TT (s) , FT (s) = TT (a) , F? (s) = J? (s) , 
FT (s) = ~\r tw - + J^T (s) , FT {s) = -^—^- r r tw - + TT {s) . (4) 
We present explicitly only one "calligraphic" quantity: 

^T(s) = ^^2 + ^MlCo(-ml-ml-s;m t ,M z ,m t )-L,(m^ 



-3B F (-s;m t ,m t )-2(l + 4r 



tz 



M 2 

2 -L ab {0,m t ,M z 



4mf — s 



(5) 



+2B F (-m 2 t ; m t , M z ) + — B F (-m 2 ; m t , M z ) + L^M 2 ) - 1 



where we introduced definitions: 

M 2 M 2 
s 2 = l _ c 2 c 2 = YL R = £ 



m 2 M 2 



/ M 2 \ 

L ab (M u M 2 , M 3 ) = ^1 + j±j M 2 C (-ml -m 2 , -s; M 3 , M 2 , M 3 ) 

-B F (-s; M 3 , M 3 ) + B F (-m 2 ; M 2 , M 3 ) , (6) 



B {-s;M 1 ,M 2 ) = - + B*{-s;M 1 ,M 2 ). (7) 

We emphasize that we leave t'Hooft scale parameter // in our formulae un- 
fixed, retaining an opportunity to control /i-independence (and therefore UV- 
finiteness) in numerical realization of one-loop form factors, providing thereby 
an additional cross-check. 

In jS] we present many examples of numerics which exhibit very good 
level of agreement between ZFITTER and new code. 



4 "Algebraic" approach to multi-loops 

There is a lot of algebraic structure in Feynman diagrams, and the idea is 
to exploit it to the maximum. That this is indeed possible was discovered 
in the so-called integration-by-parts, i-b-p, algorithm EH|- ^ P rove d to 
be hugely successful — many well-known NNLO calculations in QCD rely 
on i-b-p. In 1996, F.Tkachov came up with a mathematical result and a 
scenario to attack arbitrary multi-loop diagrams ^3]. It is far from being 
clear whether it is possible to use it even for simplest 2-loop diagrams. First 
of all, worth trying it for a familiar one-loop setting. We did it for the scalar 
form factors ^F2% iD {s) of Z cluster in the process e + e~ — ► ti, see Fig. ^ Here 
are some results we have recently obtained |15j . 

4.1 Passarino— Veltman Cq in variables of simplex 

We begin with the usual PV function Co, see, for instance [3]: 



in 2 C (pl,pl,Q 2 ;m 1 ,m,m 2 ) = fi 4 n [ d n q—\- 

J «o«i 



i 



where 



d = q 2 + ml - ie , d 1 = (q + Pif + m 2 - ie , d 2 = (q + Q) 2 + m\- ie, 
Q = Pi+P2, Q 2 = (P1+P2) 2 = s. (9) 

It is very useful to write it down in the so-called simplex variables. After 
standard calculations, we arrive at an integral over two Feynman parameters: 

C (p 2 1 ,p 2 2 ,Q 2 ;m 1 ,m,m 2 ) = I dxf dy — !— , (10) 

J0 J0 * x,y 
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with 



o 2 2 2 2 i I r^2 2 2\ , ( 2 , 2 2\ , / 2 , 2 2\ , 2 

Pxy = —x P\—y p 2 +xy(Q —p 1 —p 2 )+x{p 1 +m 1 —m ) +y [p 2 +m 2 —m )+m . 

We emphasize that both the integration domain (unit triangle) and integrand 
are symmetric with respect to interchange of x «-> y (and pi ^ p 2 , tui «-> 
m 2 ). 

Imagine now that the integral of Eq. (jl()|) is incalculable analytically (di- 
logariphms are not discovered!) and that the only way to compute it is a 
numerical integration over Feynman parameters. Representations of Eq. fjlOj) 
is very bad for numerical treatment due to existence of zeroes in polynomial 
P 



4.2 "Lifting" of polynomial powers 



As was proved in |14j . there exists a differential operator with an aid of which 
it is always possible to "lift" the power of a polynomial: 



x.y 



1 

A 



jx + A x ) d x + (y + A y ) d y 
2(k + l) 



Dk+l 
x,y ■ 



(12) 



Repeating procedure recursively several times, it is possible to transform 
any negative power into any positive power. Eq. (jl2j) involves some new 
determinant A. For instance, for P x ^ y one has: 



Px, y = Px, y (m t ,M z ,m t ) 
P X:V (m t , M z , m t ) 



Q 2 xy + m\ (x + yf + M\ (1 - x - y) 
A(m t ,M z ,m t ) 



A, 



A 



Q 



2 Ml 



4A, 



(13) 



where 

A3 = — \Q 2 {Q 2 + 4m 2 ) is the Gram determinant for our 3-point function. 
Vector A and new determinant A are inherent to a given diagram, contrary 
to Gram determinant which depends only on external momenta, they feel all 
internal masses of a diagram. 



4.3 New reduction 

Exploiting identities Eq. (jl2J) and a reduction in n-dimensions which is an 
alternative to PV reduction, we express the scalar form factors of our Z 
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cluster in terms of untaken integrals over Feynman parameters. The latter 
appear inside and outside polynomials (in powers k + e, k = 1,0). All x 
and y, which appear outside, may be eliminated yet in n-dimensions. This 
is a reduction because it reduces expressions to a very limited number of 
functions which might be termed as new scalars. Reduction in n-dimensions 
heavily exploits i-b-p and the symmetry of simplex. 

4.3.1 New scalars 

In framework of this approach we meet an analog of the usual Co and Bq 
functions: 

/l rl-x p 

dx J dyP k y \n^f, (14) 

B {k,fi 2 ;Q 2 ;M 1 ,M 2 ) = [ dxP*ln^ , (15) 

Jo ^ 

where in the most general case: P x = Q 2 x (1 — x) + M\x + M| (1 — x) . 

All new scalars, not only B , depend on the t'Hooft scale \i. We use short 
hand notation for new scalars: 



L k (m t , M z , m t ) = M z 2k C (k, fi 2 ; -m 2 , -m 2 , -s; m t , M z , m t ) , 
L k (m t ,0,m t ) = M~ 2k B (k,n 2 ;-s;m t ,m t ) , 

L k (m t ,M z ,0) = M; 2k Bo(k,fx 2 ;-m 2 ;m t ,M z ), (16) 



Next, we present some results which were derived after application of the 
procedure one and two times. 

4.3.2 "Once lifted" expressions 

We show only one scalar form factor Fl\ as a typical example: 



Fix = ^<2Lo(m t} M z ,m t )-L ll (m t ) + l L (m 4 ,M z ,0) 

c 2 w { r tz l 

-L M (M z ) + l] - ^- (L (m t ,0,m 4 ) - 2L Q (m t , M z ,m t ) - l)J. (17) 

Expression for F^f is remarkably compact, cf. Eq.fjSJ), and does not contain 
any determinants at all (the latter property is not typical, however). 
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"Twice lifted" expressions have a similar structure, however, expres- 
sions blow up with increasing the number of recursions. 

In ^Hj we show results of numerical computation of F2%, i.e. standard 
approach giving an analytic expressions Eq.fJHJ) in terms of dilogs; and of 
numerical computation of F^f, i.e. "once lifted" expressions Eq.(fT7|) and 
F2%, i.e. "twice lifted". We also show numbers for Fq Z {S12) . Numbers show 
9-10 digit agreement for Re part and 5 digits — for Im part at k — 1. This 
demonstrates that approach succeeds at one loop. 

5 Outlook 

In my opinion, nowadays in the field of theoretical support of experiments, 
it is reasonable to work in two directions: 

1) Complete one- loop SM corrections for processes with number of par- 
ticles in the final state Nf > 3 (for applications at TEVATRON, LHC, LC, 
//-factory). Our project CalcPHEP belongs to this direction. In connection 
with this, one should say that if one wants to attack a problem of a certain 
level complexity, not starting just from the previous level, all the interme- 
diate levels have to be worked through. This will considerably delay the 
realization of the project. Furthermore, complete one-loop corrections for 
the process e + e _ — > ffj (where / is any fermion, including electron) is a 
part of the two- loop program for the Z resonance, so we naturally come to 
the second direction, namely: 

2) Work towards two- loop precision level control of HEP observables (for 
applications at GigaZ LC option). "Algebraic" approach for two- loops be- 
longs to this direction. It undoubtedly works and possesses appealing features 
at one loop. Beyond one loop, there are huge algebraic difficulties; explicit 
solution has never been found. F.Tkachov currently explores scenarios that 
look bizarre enough to hold promises for such a complicated problem. Are 
there any chances of success in a foreseeable future? One may count only 
on: 1) a non-standard F.Tkachov's expertise which integrates a non-trivial 
understanding of the math involved with a considerable skill in algorithm 
design and software engineering; and on 2) our community intentions to con- 
tinue digging in this direction, in particular, now in the two-loop field (see 
Giampiero Passarino's talk at this Symposium). 

Happy birthday, Alberto! 
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